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Abstract 

Background: Polyphenol oxidase (PPO), often encoded by a multi-gene family, causes oxidative browning, a significant 
problem in many food products. Low-browning potatoes were produced previously through suppression of PPO gene 
expression, but the contribution of individual PPO gene isoform to the oxidative browning process was unknown. Here 
we investigated the contributions of different PPO genes to total PPO protein activity, and the correlations between PPO 
protein level, PPO activity and tuber tissue browning potential by suppression of all previously characterized potato PPO 
genes, both individually and in combination using artificial microRNAs (amiRNAs) technology. 

Results: Survey of the potato genome database revealed 9 PPO-like gene models, named StuPPOl to StuPP09 in this 
report. StuPPOl, StuPP02, StuPP03 and StuPP04 are allelic to the characterized P0TP1/P2, POT32, POT33 and POT72, 
respectively. Fewer ESTs were found to support the transcriptions of StuPP05 to StuPP08. StuPP09 related ESTs 
were expressed at significant higher levels in pathogen-infected potato tissues. A series of browning phenotypes 
were obtained by suppressing StuPPOl to StuPP04 genes alone and in combination. Down-regulation of one or 
several of the PPO genes did not usually cause up-regulation of the other PPO genes in the transgenic potato 
tubers, but resulted in reduced PPO protein levels. The different PPO genes did not contribute equally to the 
total PPO protein content in the tuber tissues, with StuPP02 accounting for ~ 55% as the major contributor, 
followed by StuPPOl, ~ 25-30% and StuPP03 and StuPP04 together with less than 15%. Strongly positive correlations 
between PPO protein level, PPO activity and browning potential were demonstrated in our analysis. Low PPO activity and 
low-browning potatoes were produced by simultaneous down-regulation of StuPP02 to StuPP04, but the greatest 
reduction occurred when StuPPOl to StuPPOA were all suppressed. 

Conclusion: StuPPOl to StuPP04 genes contributed to browning reactions in tuber tissues but their effect was 
not equal. Different PPO genes may be regulated independently reflecting their diversified functions. Our results 
show that amiRNAs can be used to suppress closely related members of highly conserved multi-gene family. 
This approach also suggests a new strategy for breeding low-browning crops using small DNA inserts. 

Keywords: Artificial microRNA (amiRNA), Enzymatic browning, Polyphenol oxidase gene family, Solanum 
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Background 

Polyphenol oxidase (PPO) is nearly ubiquitous in angio- 
sperms and belongs to a class of copper-binding enzymes 
that catalyze the oxidation of phenolics to quinones. The 
subsequent non-enzymatic polymerization of the qui- 
nones leads to formation of brown pigments that are 
the cause of post-harvest deterioration and loss of 
quality in many economically important crops [1,2]. 
Losses caused by the browning resulting from PPO 
catalyzed-oxidations probably account for 50% of the 
losses of industrial production of fruits and vegetables 
[3]. Several reports have described reduced browning 
reaction in crops by suppression of PPO gene expres- 
sion using transgenic transformation with PPO gene 
fragments in configurations such as sense, antisense or 
double-stranded RNA [4-10]. Those approaches func- 
tioned by establishing an RNA silencing mechanism 
guided by a population of heterogeneous small inter- 
fering RNAs (siRNAs) [11]. Inevitably, the whole PPO 
gene family in the transgenic hosts was targeted be- 
cause PPO genes are members of a multi-gene family 
with highly conserved gene sequences [2]. Because of 
this, it has been difficult to assess the contribution 
made by the individual PPO gene(s) to the oxidative 
browning process in different tissues. In addition to 
the undesired browning activity, PPOs appear to play 
important roles in signal transduction, stress and defense 
response throughout plant growth and development, but 
the specific PPO gene members involved in the different 
functions has not been elucidated. In potato (Solarium 
tuberosum L.), five PPO genes, namely POTP1 (GenBank 
ID: M95196), POTP2 (M95197), POT32 (U22921), POT33 
(U22922) and POT72 (U22923), were previously identified 
and characterized [1,12]. The nucleotide sequences of 
POTP1 and POTP2 are over 97% identical. POTP1/P2, 
POT32, POT33 and POT72 share 70-82% nucleotide hom- 
ology. A previous study based on RNA Northern blot 
analysis revealed that POTP1 and POTP2 genes were 
expressed mainly in potato leaves and flowers, POT32 
and POT33 mRNA were detected mainly in tubers 
with the POT32 being the major form throughout 
tuber development, and POT72 gene was mainly expressed 
in the roots [1]. 

Artificial microRNA (amiRNA) technology is a newly 
developed approach for inducing loss of gene function 
in plants [13-16]. It utilizes microRNA (miRNA) gene 
backbones to express artificial small RNAs that are usu- 
ally 21 nucleotides (nt) in length. The resultant amiR- 
NAs join in RNA silencing pathways and guide silencing 
of the gene of interest [17]. One of the advantages of 
amiRNA strategy is that it generates a single type of 
small RNA population all with the same selective nucleic 
acid sequence. It provides a feasible method for either si- 
lencing an individual gene or simultaneously silencing or 



partially silencing a multi-gene family while at the same 
time minimizing the risk of unpredicted off-target effects 
[18]. The amiRNA strategy has been applied in functional 
genetics studies using model plants, such as Arabidopsis 
and also agricultural crops, such as rice, alfalfa and tomato 
in recent years [19-21]. However, there were fewer reports 
of the targeting of closely related members of multi-gene 
families [20,22]. 

Here we reported suppression of the characterized mem- 
bers of the PPO gene family, i.e. POTP1/P2, POT32, 
POT33 and POT72, individually or in combination in po- 
tato using amiRNAs. This allowed us to investigate the 
contributions of the different PPO genes to the total PPO 
protein activity in potato tubers and to further understand 
the correlations between PPO protein level, PPO activity 
and tuber tissue browning potential. Our results show that 
amiRNAs can be applied to suppress the expression of in- 
dividual members of a highly conserved gene family. A 
series of browning phenotypes resulted from the suppres- 
sion of different PPO gene isoforms in potato. Our results 
also suggest a new strategy for developing low-browning 
or non-browning crops. 

The PPO gene suppression research in this report 
started before the availability of the potato genome se- 
quence data, but we recently surveyed the PPO gene 
family in the potato genome and discovered 9 PPO-like 
gene models. The PPO gene models are systematically 
named as StuPPOl to StuPP09 (Solarium tuberosum 
polyphenol oxidase 1 to 9). POTP1/P2, POT32, POT33 
and POT72 are considered allelic to StuPPOl, StuPP02, 
StuPP03 and StuPP04, respectively. For continuity in 
the systematic nomenclature in this report, we renamed 
POTP1/P2, POT32, POT33 and POT72 to StuPPOl, 
StuPP02, StuPP03 and StuPP04. 

Results 

Genome-wide survey of PPO-like gene models in 
Solanum tuberosum 

A genome-wide search of the recent S. tuberosum whole 
genome database in the US Joint Genome Institute 
(http://www.jgi.doe.gov) reveals 9 PPO-like gene models, 
tentatively named StuPPOl to StuPP09 in this report 
(Table 1, Additional file 1). The StuPPOl to StuPPOS 
genes are aligned on chromosome 8, and StuPP09 is lo- 
cated on chromosome 2. StuPPOl and StuPP06 are in 
close proximity to each other in a 47-kb region, while 
StuPP02, StuPP03, StuPP04, StuPPOS, StuPP07 and 
StuPP08 are clustered in tandem in a 144-kb region. 
The two regions are separated by a distance of 2,072-kb 
on chromosome 8. Analysis of the deduced amino acid 
sequences of the major peptides encoded by the PPO-like 
genes suggests that the predicted proteins all contain three 
typical PPO protein domains: the tyrosinase domain 
(pfam00264), the PPO!_DWL domain (pfaml2142) and 
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Table 1 List of predicted potato PPO gene models 


Tentative gene 
name (in this report) 


Locus location 


Number of 
intron 


Predicted transcript name 
assigned by potato genome 
sequencing consortium (PSGC) 


Possible Allele (GenBank ID, 
nucleotide sequence identity 0 /)) 


StuPPOl 


chr08, 30458794.30460741 


0 


PSGC003DMT400076054 


P0TP1 (M95196, 97 .3%)/ 'P0TP2 (M95197, 

r\~i nn/ \ /vn i\ r\r\r o r r 1 t'\r\r\c\/\ 

y/.y%j/x(vi_00D3bb i // ( i oo%) 


C+, , DD/O") 

btuPPUz 


cnruo, ozo/z i y4..ozo/4 1 yz 


0 


rbuLUUUoDIVI 1 4UUU4ooo4 


rUlDZ (UZZyzl, yO./%J/XIVI_UUOOOJOZ 1 (yj.JTOj 


JLUr r KJj 


LMlUO, jZDO/ JjU..jZOoyZOU 


n 
u 




\\J\JJ yJZZyZZ, yH./ /Oj/ AIVI_UUDjDJ jZU 

(100%) 


StuPP04 


chr08, 32667904..32669792 


0 


PSGC0003DMT400048685 


P0772 (U22923, 96.8%) 


StuPP05 


chr08, 32591 830..32593339 


0 


PSGC0003DMT400048692 


XR_1 83056 (84.6%) 


StuPP06 


chr08, 30504049..30505788 


0 


PSGC0003DMT400076055 


XM_004245989 (89.1%) 


SWPP07 


chr08, 32577708..32579291 


0 


PSGC0003DMT400048703 


XR_1 83056 (84.8%) 


StuPP08 


chr08, 3270381 8..32721 729 


2 


PGSC0003DMT400048679 


XM_006365329 (100%) 


StuPP09 


chr02, 5559371 8..5559601 9 


1 


not available 


XM_006347021 (100%) 



the PP01_KFDV domain (pfaml2143) [23,24], but the 
tyrosinase domains with the putative StuPP05 and 
StuPP07 peptides are shorter than the others and in- 
complete (Additional file 1: Part B). StuPPOl, StuPP02, 
StuPP03 and StuPP04 are possibly allelic to the charac- 
terized POTP1/P2, POT32, POT33 and POT72, re- 
spectively, considering that the nucleotide sequences 
between the potential alleles are 95-99% identical 
(Table 1, Additional file 1). Numerous ESTs were 
found from different developmental potato tissues 
for StuPPOl to StuPP04 (Additional file 1: Part A 
and Part C). Surprisingly, StuPPOl appears to be the 
only possible allele to the POTP1 and POTP2 genes, 
and no duplication of the StuPPOl locus was ob- 
served by searching the S. tuberosum genome se- 
quence. StuPPOS, StuPP06 and StuPP07 are three 
novel PPO-like gene models predicted from this S. 
tuberosum genome sequence analysis. However, only 
few EST fragments (0 to 3) that probably relate to the 
potential transcripts of these three genes were found, and 
the ESTs cover only fragmental regions of the putative 
transcripts (Additional file 1: Part A and Part C). The 
low prevalence in EST databases suggests that StuPPOS, 
StuPP06 and StuPP07 may be expressed at very low levels 
in S tuberosum. StuPPOS and StuPP09 are the only PPO- 
like gene models with introns. No EST from S. tuberosum 
EST databases was found for StuPP08, suggesting that this 
gene sequence is normally not transcribed. StuPP09 is the 
only PPO-like gene model that is not clustered with the 
others on chromosome 8, but is independently located on 
chromosome 2. A number of ESTs were found for 
StuPP09, but all of the ESTs were revealed in the cDNA 
libraries from the tissues of in vitro cultured potato callus 
(DBLINK ID: LIBEST_0 15047), abiotic stress treated leaf 
and root (LIBEST_015048) [25], or pathogen-infected leaf 
and tuber (LIBEST_008810, LIBEST_009854, LIBEST_ 
015324, LIBEST_015920, LIBEST_017649, and LIBEST_ 



025550) (Additional file 1: Part A and Part C) [26-28]. The 
expression data of the supporting ESTs of StuPP09 were 
mostly not available because the gene model was not anno- 
tated previously and most reports focused on annotated 
genes. However, at least three of the StuPP09 related ESTs 
(GenBank ID: CK640809, CO267905 and GT888802) were 
found to express differentially in the pathogen-infected po- 
tato tissues (Additional file 1: Part C). For example, the ex- 
pression level of the CK640809 was induced 3- to 14-fold 
higher in potato (cultivar Indira and Bettina) leaf tissue that 
was inoculated with fungus Phytophthora infestans (See the 
Table four in reference [27]. The CO267905 showed over 
2-fold induction in the potato (cv. R-gene-free potato clone 
386209.10) leaf tissue infected with P. infestans at 24 hours 
post-inoculation (hpi) and its expression level was over 17- 
fold higher at 48 hpi (See the Table one in reference [28]). 
The relative expression level of GT888802 was about 3-fold 
higher in potato (cv. Spunta) tubers inoculated with fun- 
gus Fusarium eumartii at 24 hpi (Table S2 in reference 
[26]). These data imply that StuPP09 is probably an in- 
ducible PPO gene expressed in response to disease 
defense and cell rescue [28]. 

Generation and selection of amiRNA-expressed transgenic 
potatoes 

Seven different amiRNAs were designed to target the 
characterized PPO genes, namely StuPPOl (previously 
named POTPl/P2\ StuPP02 (previously, POT32), StuPP03 
(previously, POT33) and StuPP04 (previously, POT72) in 
S. tuberosum [1,12], with amiRNA sequences comple- 
mentary to either a specific gene or multiple targets by 
choosing the appropriate 21 -bp region of the corre- 
sponding PPO genes (Table 2). The amiRNA sequences 
were incorporated in an Arabidopsis thaliana miR168a 
gene backbone built in a plant binary vector (Figure 1), 
[29,30]. From transformation of thousands of explants, 
8 to 10 transgenic potato lines for each amiRNA 
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Table 2 List of the amiRNA constructs, amiRNA names, amiRNA target sequences and target genes 



Construct 


amiRNA name 


amiRNA sequence (5'— > 3') 


Target sequence* (5'— > 3') 


Gene name (GenBank accession ID) 


prZamiKrrU I 


aminrrU 1 


i ii \r~(~\ \r i\r\ \rr\ \rri\i\\ n \c at 


c\ \c a ai ii \cc hrc An irArrAA 
UULAAUUULALLAUULALLAA 


JLUrrU 1/rUI r / IrZ (AIVI_UUodjj 1 /// 










M95196/M95197) 


pPZamiRPP02 


amiRPPP02 


U UGCU AGCUGGCGGAAGUGAA 


UUCACUUCCGCCAGCUAGCAA 


StuPP02/POT32 (U22921) 


pPZamiRPP03 


amiRPP03 


UUGUUCACUGGGGGGAGUGUA 


UACACUCCCCCCAGUGAACAA 


POT33 (U22922) 








UUCACUCCCCCCAGUGAACGA 


StuPP03 (XM_006365320) 


pPZamiRPP023 


amiRPP023 


UCAUCAACUGGAGUUGAGUUG 


CAACUCAACUCCAGUUGAUGA 


StuPP02/POT32 








CAACUCAACUCCAGUUGAUGA 


StuPP03/POT33 


pPZamiRPP0234 


amiRPP0234 


UAGAACUCGGAGUUCAACCAA 


UUGGUUGAACUCCGAGUUCUA 










UUGGUUGAACUCCGAGUUCUU 


SWPP02/POT32 








UUGGUUGAACUCCGAGUUCUU 


StuPP03 (XM_006365320) 








UUGGUUGAACUCUGAGUUCUU 


POT33 (U22922) 








UUGGUUGAACUCCGAGUUCUU 


SWPP04/POT72 (U 22923) 


pPZamiRPP0234A 


amiRPP0234A 


AAGAACUCGGAGU UCAACCAA 


UUGGUUGAACUCCGAGUUCUU 


StuPP02/POT32 








UUGGUUGAACUCUGAGUUCUU 


SWPP03/POT33 








UUGGUUGAACUCCGAGUUCUU 


SWPP04/POT72 


pPZamiRPP01234 


amiRPP01234 


UCAAGCUCAUUCGCAUUCACA 


UGUGAAUGCGAAUGAGCUUGA 










UGUGAAUGCGGAUGAGCUUGA 


StuPPO 1/P0TP 1/P2 








UGUGAAUGCAAAUGAGCUUGA 


StuPP02 (XM_006365321) 








UGUGAAUGCGAAUGAGCUUGA 


POT32 (U22921) 








UGUGAAUGCGAAUGAGCUUGA 


SWPP03/POT33 








UGUGAAUGCGAAUGAGCUUGA 


StuPP04/POT72 



Note: * The non-identical nucleotides, compared to the target sequence, are marked in bold. The homologous sequences of the predicated PPO-like gene models 
{StuPPO! to SWPP09) to all amiRNA target sequences are aligned in the Additional file 1: Part D. 



miR16Sa 



5" jTCGCTTGGTGCAGGTCGGGAA | [ ATCCCGCCTTGCATCAACTGA 



amiRPPOl 




TTGGTGACTGGTGCAATTGAC 



amiRPPOr 
] [ GTCAATTGCGCCAGTCTCCAT 



amiRPPQ2 amiRPP02* 

TTGCTAGCTGGCGGAAGTGAA | [ TT C ACTT C C A C C A G CTT G CAT 



amiRPP03 



amiRPPOS* 



TTGTTCACTGGGGGGAGTGTA | 1 TA C A CTC C CTC C AGT G TAC AT 

amiRPP023 amiRPPQ23« 

] 1 CAACTCAACGCCAGTTAATGT 



TCATCAACTGGAGTTGAGTTG 



amiRPP0234 



amiRPPQ234* 



TAG A AC T C G G A G TTC A AC C A A | [ TTGGTTGAATTCCGAGCTCTT 



amiRPP0234A amiRPPQ234A" 

A AG A ACTC G G A G T T C A ACC AA | [ T TG GTT G A ATT C C G AGC TCTA 



amiRPP01234 



TCAAGCTCATTCGC ATTCACA 



amiRPPOl 234* 

] [ TGTG AATG C A A ATG AGTTTGT 



- |35S-T| // 

H 36S T | // 



MiR ;68a(EU549054. 1) 

B 

pPZamiRPPOl 
pPZamiRPP02 
pPZarniRPPQ3 
pPZamiRPP023 
pPZamiRPP0234 
P PZamiRPP0234A 
pPZamiRPP01234 

Figure 1 Diagrammatic representation of artificial microRNA constructs. (A) Linear structure of the miR168a primary transcript gene 
(MIR 168a, nt 120 to 355, GenBank Accession No. EU549054.1). Sequences of the miR168a and its complementary region (illustrated as 
approximately miR168a*) in the gene are displayed in the boxes. (B) Structure of the binary vectors for expression of amiRNAs. Construct names 
are indicated at the left. The sequences of the designed amiRNA and its complementary region (approximately amiRNA*) are displayed in the 
boxes. 35S-P, CaMV 35S promoter. 35S-T, CaMV 35S terminator. 



■3' 



Hass-Tl — // 

H 35S T [ // 



{35S 



5W/ 
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construct were selected and propagated for molecular 
genetic screening and analysis (Table 3). Significant 
down-regulation of the targeted gene expression was 
detected in a number of the resulting transgenic potato 
lines (Table 3). Based on the initial real-time quantita- 
tive reverse transcription PCR (qRT-PCR) assay of the 
in vitro cultured potato plants, the transcript levels of 
the StuPPOl genes were reduced by 68 to 98% in six 
amiRPPOl series transgenic lines (clones). The amiRPPOl 
series expressed artificial miRNA - amiRPPOl designed to 
target the StuPPOl gene. Similarly, one amiRPP02 series 
line (target: StuPP02), two amiRPP03 series lines (target: 
StuPP03), five amiRPP023 series lines (targets: StuPP02 
and StuPP03), two amiRPP0234 series lines (targets: 
StuPP02, StuPP03 and StuPP04\ three amiRPP0234A 
series lines (targets: StuPP02, StuPP03 and StuPP04) and 
five amiRPP01234 series lines (targets: StuPPOl, StuPP02, 
StuPP03 and StuPP04) showed substantial reduction of 
the target gene transcript level(s) (Table 3). Mature amiR- 
NAs were detected by reverse transcription PCR (RT- 
PCR) in the selected amiRPPOl, amiRPP03, amiRPP023, 
amiRPP0234A and amiRPP01234 series lines (Table 3, 
Additional file 2: Figure SI). However, no mature amiR- 
NAs were revealed by RT-PCR in the 10 lines of the 
amiRPP02 series nor the 10 clones of the amiRPP0234 
series (Table 3 and data not shown). The amiRPP02- and 
amiRPP0234- transgenic lines were therefore excluded 
from further evaluations. In addition, small RNA North- 
ern blots were previously done to detect amiRNA expres- 
sion of the multiple transgenic lines listed in Table 3, 
including amiRPPOl -7 and -12 (previously named as 
amiR-POTPl/P2, L7 and L12), amiRPP02-12, -19 (previ- 
ously, amiR-POT32, L12 and L19), amiRPP03-8 and -15 
(previously, amiR-POT33, L8 and L15) and amiRPP0234- 
9 and -10 (previously, amiR-POT32/33/72, L9 and L10) [30]. 

In plants, miRNA-guided RNA silencing has been shown 
to occur mostly through complementary cleavage of the 
targeted mRNA by Argonaute proteins [31]. Considering 
the highly complementary sequences between the designed 
amiRNAs and their target PPO genes, we used 5 '-RACE 
PCR to detect the possible cleavage of the PPO gene tran- 
scripts). Because of the multiple PPO gene members in 
potatoes, we developed a strategy for detecting all pos- 
sible fragmented-mRNA of the PPO genes but not the 
5 '-capped mRNA (Figure 2A and see the Methods). As 
predicted, a 253-bp PCR product (including a 45-bp 
5 '-RACE adaptor) was revealed by the nested PCR 
round- 1 from the enriched poly(A) + RNA of the young 
leaves of line amiRPPOl-12 (Figure 2C). Sequence analysis 
demonstrated the fragment included two nearly identical 
sequences differing by one nucleotide (Additional file 3: 
Figure S2). Both were highly related to the StuPPOl gene 
(97 to 99% identity at nucleotide level) (Additional file 4: 
Figure S3). The first 10 nucleotides of the 5 '-end of the 



sequences were complementary to the 5 '-end of the 
designed amiRPPOl, indicating the products were 
from cleavage of the target StuPPOl mRNA and the 
cleavage site was between nucleotides 10 and 11 at the 
amiRPPOl site (Figure 2E). The presence of the cleaved 
StuPPOl mRNA was also demonstrated by the specific 
nested PCR-2 (Figure 2D). The results suggested that the 
expressed amiRNAs in the transgenic plants functioned as 
the small RNAs that determined the silencing of the gene 
(s) of interest. The following transgenic lines were selected 
and propagated for further biological analysis, amiRPPOl - 
2, -3 and -12, amiRPP03-12 and -15, amiRPP023-5, -7 
and -9, amiRPP0234A-4, -6 and -14, amiRPP01234- 
2, -6 and -12 (Table 3). 

PPO gene expressions in transgenic potatoes 

No growth abnormities occurred in the amiRNA-expressed 
transgenic plants under greenhouse conditions. Nor were 
the range of the tuber sizes and weights significantly differ- 
ent in either the transgenic or the wild types (Data not 
shown). Relative transcript levels of StuPPOl, StuPP02, 
StuPP03 and StuPP04 genes in tuber tissues of the trans- 
genic lines were assayed by qRT-PCR with the results illus- 
trated in Figure 3. For lines amiRPPOl -2, -3 and -12, the 
expression of the target StuPPOl gene was suppressed by 
90 to 99%. The mRNA levels of the non-targeted StuPP02 
and StuPP03 in lines amiRPPOl-3 and amiRPPOl-12 
were similar to that of the non-transgenic, wild types 
WTJ potato controls, but the expression levels of the two 
genes (StuPP02 and StuPP03) in line amiRPPOl-2 
were unexpectedly reduced by 50% and 80%, respect- 
ively (Figure 3A). For lines amiRPP03-12 and -15, the 
transcript of the target StuPP03 gene was reduced by 
over 75%. In addition, the non-targeted StuPPOl and 
StuPP02 gene mRNA levels also decreased by 50-60% in 
line amiRPP03-12, but the two gene transcripts {StuPPOl 
and StuPP02) in line amiRPP03-15 were close to the level 
observed in the WT. The non-targeted StuPP04 gene 
mRNA was reduced by -60% in line amiRPP03-15 but 
the same gene mRNA level in line amiRPP03-12 was 
almost the same as the WT (Figure 3B). For lines 
amiRPP023-5, -7 and -9, the two targeted genes, StuPP02 
and StuPP03 were almost completely silenced (> 95%) and 
the non-targeted StuPPOl gene was expressed at a level 
similar to the WT, but the StuPP04 gene mRNA level 
was also generally 50-80% lower (Figure 3C). For lines 
amiRPP0234A-4, -6 and -14, the mRNA levels of all 
three targets, StuPP02, StuPP03 and StuPP04 were 
75-95% lower than the WT, but the expression of the 
non-targeted StuPPOl gene was also reduced by an 
average of 50% (Figure 3D). For lines amiRPPO 1234-2, -6 
and -12, the mRNA levels of the targeted StuPPOl, 
StuPP02, StuPP03 and StuPP04 genes were all sup- 
pressed by 85-99% (Figure 3E). 
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Table 3 Screening amiRNA-expressed transgenic potato lines 



Relative expression level of target gene 
Transgenic line Kan Insert amiRNA 

P0TP1/P2 POT32 POT33 POT72 



amiRPP01-2 


+ 


+ 


0.02 


±0 






+ 


amiRPP01-3 


+ 


+ 


0.06 


±0 






+ 


amiRPP01-7 


+ 


+ 


0.14 


±0 






+ 


amiRPP01-8 


+ 


+ 


0.03 


±0 






n.t. 


amiRPP01-9 


+ 


+ 


0.15 


±0.01 






n.t. 


amiRPPOI-10 


+ 


n.t. 


n.t. 








n.t. 


amiRPP01-12 


+ 


+ 


0.34 


±0 






+ 


amiRPP01-13 


+ 


n.t. 


n.t. 








n.t. 


WT 


_ 


_ 


0.99 


±0.08 






_ 


amiRPP02-1 


+ 


n.t. 






n.t. 




_ 


amiRPP02-4 


+ 


n.t. 






1.42 ±0.12 




_ 


amiRPP02-9 


+ 


+ 






n.t. 




_ 


amiRPP02-12 


+ 


+ 






1 .03 ± 0.02 




_ 


amiRPP02-13 


+ 


n.t. 






n.t. 




_ 


amiRPP02-15 


+ 


n.t. 






0.83 ±0.11 




_ 


amiRPP02-16 


+ 


+ 






0.45 ± 0.03 




_ 


amiRPP02-17 


+ 


n.t. 






n.t. 




_ 


amiRPP02-18 


+ 


+ 






1 .07 ± 0.04 




_ 


amiRPP02-19 


+ 


+ 






0.99 ± 0.03 




_ 


WT 


_ 


_ 






1.00 ±0.01 




_ 


amiRPP03-8 


+ 


+ 








1 .29 ± 0.23 


+ 


amiRPP03-9 


+ 


n.t. 








n.t. 


n.t. 


amiRPP03-12 


+ 


+ 








0.66 ±0.1 2 


+ 


amiRPP03-13 


+ 


n.t. 








0.97 ±0.1 7 


n.t. 


amiRPP03-14 


+ 


+ 








n.t. 


n.t. 


amiRPP03-15 


+ 


+ 








n.v. 


+ 


amiRPP03-16 


+ 


n.t. 








1.31 ±0.40 


n.t. 


amiRPP03-17 


+ 


+ 








n.t. 


n.t. 


amiRPP03-18 


+ 


+ 








n.t. 


n.t. 


WT 


_ 


_ 








1 .00 ± 0 


_ 


amiRPP023-3 


+ 


+ 






n.t. 


n.t. 


n.t. 


amiRPP023-4 


+ 


n.t. 






1.03 ±0.11 


1 .28 ± 0.29 


n.t. 


amiRPP023-5 


+ 


+ 






0.32 ±0.01 


0.44 ± 0.08 


+ 


amiRPP023-7 


+ 


+ 






0.28 ± 0.03 


0.19 ± 0.01 




amiRPP023-8 


+ 


+ 






0.42 ± 0 


0.01 ±0 


+ 


amiRPP023-9 


+ 


+ 






0.07 ±0.01 


0.01 ±0 


+ 


amiRPP023-14 


+ 


n.t. 






n.t. 


n.t. 


n.t. 


amiRPP023-15 


+ 


n.t. 






0.23 ± 0.07 


0.00 ± 0 


n.t. 


amiRPP023-16 


+ 


n.t. 






n.t. 


n.t. 


n.t. 


WT 










1 .00 ± 0.08 


1.00 ±0.12 




amiRPP0234-3 


+ 


n.t. 






n.t. 


n.t. 


n.t. 


amiRPP0234-4 


+ 


n.t. 






n.t. 


n.t. 


n.t. 


amiRPP0234-5 


+ 


+ 






0.87 ± 0.06 


0.92 ±0.14 


1.03 ±0.17 
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Table 3 Screening amiRNA-expressed transgenic potato lines (Continued) 



amiRPP0234-6 


+ 


n.t. 




0.98 ±0.1 2 


3.39 ± 0.47 


1 .50 ± 0.08 




amiRPP0234-7 


+ 


+ 




1.22 ±0.16 


1.16 ±0.09 


1.04 ±0.17 




amiRPP0234-8 


+ 


+ 




1 .66 ± 0.02 


0.71 ±0.04 


0.69 ± 0.02 




amiRPP0234-9 


+ 


+ 




0.47 ± 0.04 


0.40 ± 0.09 


0.50 ±0.1 3 




amiRPPO234-10 


+ 


+ 




0.61 ±0.03 


0.61 ±0.09 


0.54 ±0.21 




amiRPP0234-1 1 


+ 


n.t. 




1.35 ±0.16 


0.83 ±0.1 7 


1.16 + 0.13 




amiRPP0234-13 


+ 


n.t. 




n.t. 


n.t. 


n.t. 




WT 








1.00 ±0.01 


1 .00 ± 0.08 


1.00 ±0.01 




amiRPP0234A-2 


+ 


n.t. 




1 .09 ± 0.20 


1.53 ±0.19 


1.26 ±0.13 


n.t. 


amiRPP0234A-4 


+ 


+ 




0.34 ±0.05 


0.43 ± 0.07 


0.04 ±0.01 


+ 


amiRPP0234A-6 


+ 


+ 




0.09 ±0.01 


0.02 ± 0 


0.05 ±0.01 


+ 


amiRPPO234A-10 


+ 


+ 




1 .07 ± 0.03 


0.77 ± 0.25 


0.47 ± 0.02 


+ 


amiRPP0234A-13 


+ 


+ 




n.t. 


n.t. 


n.t. 


n.t. 


amiRPP0234A-14 


+ 


+ 




0.58 ±0 


0.22 ±0.01 


0.1 5 ±0.03 


+ 


amiRPP0234A-15 


+ 


n.t. 




1.25 ±0.18 


0.8 ±0.1 2 


1.34 ±0.16 


n.t. 


amiRPP0234A-16 


+ 


n.t. 




0.83 ±0.1 2 


1.21 ±0.09 


0.60 ± 0.07 


n.t. 


WT 








0.99 ± 0.08 


1.01 ±0.18 


0.99 ± 0.05 




amiRPPOl 234-1 


+ 


n.t. 


n.t. 


n.t. 


n.t. 


n.t. 


n.t. 


amiRPPOl 234-2 


+ 


+ 


0.52 ± 0.07 


0.30 ± 0.03 


0.46 ±0.01 


0.34 ± 0.08 


+ 


amiRPPOl 234-3 


+ 


+ 


0.43 ± 0.04 


0.1 9 ±0.06 


0.07 ± 0 


0.43 ±0.1 


n.t. 


amiRPPOl 234-6 


+ 


+ 


0.27 ± 0.03 


0.20 ± 0.04 


0.19 + 0 


0.2 ± 0.04 


+ 


amiRPPOl 234-7 


+ 


n.t. 


n.t. 


n.t. 


n.t. 


n.t. 


n.t. 


amiRPPOl 234-8 


+ 


n.t. 


n.t. 


n.t. 


n.t. 


n.t. 


n.t. 


amiRPPOl 234-9 


+ 


+ 


n.t. 


n.t. 


n.t. 


n.t. 


n.t. 


amiRPPOl 234-11 


+ 


n.t. 


0.1 5 ±0.01 


0.34 ±0.11 


1.85 ±0.17 


0.08 ±0.01 


n.t. 


amiRPPOl 234-1 2 


+ 


+ 


0.28 ± 0.03 


0.66 ± 0.02 


0.47 ± 0.03 


0.1 8 ±0.01 


+ 


amiRPPOl 234-1 6 


+ 


n.t. 


0.42 ± 0.03 


0.1 4 ±0.03 


0.15 + 0 


0.08 ± 0.06 


n.t. 


WT 






1 .00 ± 0.09 


1.00 ±0.1 


1.00 ±0.19 


1 .00 ± 0.07 





Note: Kan, is selection of transgenic lines with Kanamycin (100 mg/L); +, Kanamycin resistant; -, Kanamycin susceptible (such as for WT). Insert, is transgene 
insertion into the chromosome(s) detected by PCR; +, positive in PCR detection; -, negative in PCR detection; n.t., not test. Relative expression level of target 
gene, was assayed using in vitro cultured potato leaf and stem tissues by qRT-PCR; the data is the average of three technical repeats and their standard deviation 
of using one original sample from each transgenic line; n.v., no value; n.t. not test. amiRNA, is detection of mature amiRNAs by RT-PCR; +, detectable (positive); 
-, not detectable (negative); n.t, not test. 



PPO protein level 

Total PPO protein levels in the transgenic and non- 
transgenic potato tubers were analyzed using a semi- 
quantitative protein dot-blot assay. Figure 4A shows the 
values of the total PPO protein level in the transgenic 
potato tuber tissues relative to that in the non-transgenic 
wild types (relative PPO protein level, simply represented 
by 'RPR' in this report). The average RPR values for lines 
amiRPPOl-2, -3 and -12 ranged from 0.70 to 0.77, indi- 
cating total PPO protein levels in these transgenic potato 
tuber tissues were on average about 23-30% lower than 
those in the wild type. Lines amiRPP03-12 and -15 
showed a decline of 15-20% on average in their total PPO 
protein level, compared to the wild type. Average reduc- 
tions of 45-70% in total PPO protein level were detected in 



lines amiRPP023-5, -7 and -9, amiRPP0234A-4, -6 and 
-14, based on their average RPRs (0.30-0.55). The average 
RPRs for lines amiRPPOl 234-2, -6 and -12 varied from 
0.20 to 0.27, suggesting that total PPO protein concentra- 
tions in these transgenic tubers decreased on average 
by 73-80%, compared to the wild type (Figure 4A). 

PPO enzymatic activity 

Figure 4B depicts the PPO enzymatic activities in the 
transgenic potato tuber tissues relative to those in the 
wild type (relative PPO activity, abbreviated to 'RPPO' in 
this report). Based on the RPPOs, the average PPO activ- 
ity of lines amiRPPOl-2, -3 and -12 was 25-35% lower 
than that of the wild type. A reduction of 15-25% in 
PPO activity was observed in lines amiRPP03-12 and -15. 
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Specific PGR primers amiRPPOl 3' CAGUUAACGUGGUCAGUGGUU 5' 



Figure 2 Detection of cleaved PPO gene mRNAs by 5' RACE-PCR. (A) Schematic diagram of a strategy for detecting the truncated but not 
the 5'-capped PPO gene mRNAs in amiRNA-expressed transgenic plants. (B) RACE-PCR (1st-round PCR). Lane 1, DNA ladder (NEB, Cat# N0474S); 
lane 2, 1st-round PCR result from the enriched young leaf poly(A) + RNA of transgenic line amiRPP01-12. The line with an arrowed end indicates 
the band of the expected size. (C) Nested PCR round-1 (using the RACE-PCR product as the template). Lane 1, DNA ladder; lane 2, nested PCR-1 
result; the line with an arrowed end indicates the band of the expected size. (D) Nested PCR round-2 (using nested PCR-1 product as template). 
Lane 1, DNA ladder; lane 2, PCR for detection of truncated StuPPOl gene mRNA; amplified band size as indicated. Lane 3, PCR for detection of 
truncated StuPP02 gene mRNA; the faint bands were non-specific amplification. Lane 4, PCR for detection of truncated StuPP03 gene mRNA; the 
faint bands were non-specific amplification. Lane 5, PCR for detection of truncated SWPP04 gene mRNA; the faint bands were non-specific ampli- 
fication. (E) Determination of the target cleavage site of amiRPPOl by sequencing the 253 bp of the nested-PCR-1 product (Figure 3C). The 
target sequences are aligned with the amiRPPOl complementarily. The arrowed line indicates that 7 of the 7 clones (7/7) were the products 
cleaved at the expected site. See Additional file 3: Figure S2 and Additional file 4: Figure S3 for the full sequence analysis. 



PPO activities in lines amiRPP023-5, -7 and -9, 
amiRPP0234A-4, -6 and -14, amiRPP01234-2, -6 
and -12 were in a similar range, about 75-95% less 
than those of the non-transgenic controls (Figure 4B). 

Browning potential and browning phenotype 

Browning was used to measure the potential of phenolic 
oxidation after mechanical release of PPO proteins from 
their storage site in plastids [32]. Browning potential of 
the transgenic potato tuber was compared to that of 
their non-transgenic wild types (relative browning po- 
tential, abbreviated to 'RBR' in this report) and the re- 
sults are shown in Figure 4C The RBRs for lines 
amiRPPOl-2, -3 and -12, amiRPP03-12 and -15 ranged 
between 0.65 and 0.90, suggesting that browning po- 
tentials of these transgenic lines were about 10-35% lower 
than those of the WT. Browning potentials of lines 



amiRPP023-5, -7 and -9, amiRPP0234A-4, -6 and -14 
were reduced by -50-65% based on their RBR values of 
~0.35-0.50. The ranges of 0.25-0.35 in RBR for lines 
amiRPP01234-2, -6 and -12 indicated that the browning 
potential of these transgenic lines was about 65-75% lower 
than the comparable wild type (Figure 4C). The browning 
potential results were relatively consistent with the series 
of visible browning phenotypes displayed after air- 
exposure of the freshly sectioned potato tubers at room 
temperature (Figure 5). Browning or blackening tissues 
developed on the sectioned tuber surfaces, typically start- 
ing from the vascular ring region and advancing to the 
medulla with increased exposure time to oxygen in the 
air. The wild type tubers developed brown tissues more 
quickly, over a larger area and more severely than in the 
transgenic tubers. Among the different transgenic types, 
the amiRPPOl and amiRPP03 series lines showed 
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Figure 3 Relative transcript levels of PPO genes in transgenic 
potato tuber tissues. (A) Transgenic lines of series amiRPPOL 
(B) Transgenic lines of series amiRPP03. (C) Transgenic lines of series 
amiRPP023. (D) Transgenic lines of series amiRPP0234A. (E) Transgenic 
lines of series amiRPP01234. Each column represents the mean value 
obtained from qRT-PCR performed in triplicate on each biological 
sample. The bars indicate standard deviation. Two biological replicates 
(indicated as a and b) from each transgenic line were selected for 
the assay. Cyclophilin and ef1a genes were used as normalization 
references and non-transgenic potatoes (WT) were set as the 
control. 



relatively stronger browning phenotype, followed by the 
amiRPP023, amiRPP0234A and amiRPP01234 series 
lines, based on the degree of dark color development that 
ranged from high to low (Figure 5). 



Statistical correlations 

Pearsons correlation coefficient (r 2 ) analysis indicated 
significantly strong and positive correlations between 
RPR, RPPO and RBR. (r 2 = 0.85-0.89, P < 0.0001) in po- 
tato tuber tissues (Table 4). Among the potato PPO 
genes, the StuPP02 gene was highly correlated with 
RPR, RPPO and RBR (r 2 = 0.70-0.80, P < 0.0001). Both 
the StuPP03 and StuPP04 genes were moderately associ- 
ated with RPR, RPPO and RBR (r 2 = 0.59-0.71, P < 0.0001), 
while the StuPPOl gene had a weak correlation with the 
three browning-related parameters, RPR, RPPO and RBR 
(r 2 = 0.19-0.27, P < 0.05) (Table 4). 

Principal component analysis (PCA) generated only 
two principle components with eigenvalues exceeding 
1.0 (Kaisers rule) (Figure 6 and Additional file 5: Table 
SI). The two components explained 87% of the total vari- 
ance. The first principle component (PCI) accounted for 
71% of total variance and had approximately equal posi- 
tive loading for the variables StuPP02 gene, StuPP03 
gene, StuPP04 gene, RPR, RPPO and RBR. Each of the 
above variables contributed about 14-18% to the PCI, sug- 
gesting their equivalent proportion in the different trans- 
genic lines. In contrast, the variable StuPPOl gene also 
contributed positively to the PCI but with a lower score 
(< 0.8% contribution). The second principle component 
(PC2) only accounted for 16% of the total variance and was 
mainly influenced by positive loading of StuPPOl gene 
(contribution to PC2, 79%) (Additional file 5: Table SI). 
The score plot of the PCI and PC2 paralleled the distribu- 
tion of the browning phenotypes (Figure 6). Lines 
amiRPPOl-2, -3 and -12, amiRPP03-12 and -15 and the 
WT, susceptible to browning were scattered on the right 
side of the plot. In contrast, lines amiRPP023-5, -7 and 
-9, amiRPP0234A-4, -6 and -14, amiRPPO 1234-2, -6 
and -12, resistant to browning were separated to the left 
side of the plot. Noticeably, lines amiRPPO 1234-2, -6 and 
-12, which inhibited StuPPOl, StuPP02 } StuPP03 and 
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Figure 4 Assay of PPO protein level, enzymatic activity and 
browning potential in transgenic potato tuber tissues. (A) Relative 
PPO protein level (RPR). (B) Relative PPO activity (RPPO). (C) Relative 
browning potential (RBR). Each box plot presents the data from three 
biological and three technical replicates of the transgenic and 
non-transgenic (WT) potato tubers. All data are presented relative 
to the level of the WT. A line across the box indicates the median. 
The box indicates the 25th and 75th percentiles. Whiskers represent the 
maximum and minimum values. Different lower case letters indicate 
values are significantly different at P < 0.05 level; different capital 
letters indicate values are highly significantly different at P < 0.01 level 
based on Duncan's Multiple Range Test. 



Hierarchical clustering analysis of the transgenic lines on 
the variables RPR, RPPO, RBR and PPO gene expression 
levels produced two major clusters, subcluster-1 (top) and 
subcluster-2 (bottom) (Figure 7). Members in each of the 
subclusters displayed a similar pattern with regards to ex- 
pression trends in the variables RPR, RPPO and RBR. Statis- 
tically, subcluster-1 expressed considerably higher values 
(Min, Max and range) of RPR, RPPO and RBR than 
subcluster-2 (Figure 4). Interestingly, subcluster-1 consisted 
the WT and the transgenic lines designed for targeting a sin- 
gle PPO gene (StuPPOl or StuPP03\ including lines 
amiRPPOl-2, -3 and -12, amiRPP03-12 and -15, whereas 
subcluster-2 included transgenic lines designed for targeting 
multiple PPO genes, namely, the lines amiRPP023-5, -7 
and -9, amiRPP0234A-4, -6 and -14, amiRPP01234- 
2, -6 and -12. The transgenic lines were further di- 
vided into smaller sub-groups based on their different 
scores (Figure 7). For example, the group of lines 
amiRPPOl-3 and amiRPPOl-12 displayed similar level 
of StuPPOl gene suppression but showed almost nor- 
mal gene transcript levels of StuPP02, StuPP03 and 
StuPP04 (Figure 3A). Although the gene transcript 
levels of both StuPP03 and StuPP04 in line 
amiRPP03-15 were lower by -70% than those of the 
WT (Figure 3B), the two clustered together because 
the values of RPR, RPPO and RBR in amiRPP03-15 
were closer to those in the WT than other transgenic 
lines (Figure 4A, B and C). Lines amiRPPOl-2 and 
amiRPP03-12 grouped based on their similar gene ex- 
pression levels of StuPP02, StuPP03 and StuPP04 
(Figure 3A and B) and the generally similar values of 
RPR, RPPO and RBR (Figure 4A, B and C). The group of 
lines amiRPP023-7 and amiRPP023-9 showed strong 
down-regulation of both the StuPP02 and StuPP03 genes, 
moderate down-regulation of the StuPP04 gene (Fig- 
ure 3C), as well as similar RPPO and RBR values (Fig- 
ure 4B and C). Three amiRPP0234A series lines (-4, 
-6 and -14) clustered together based on their similar 
down-regulated gene expression levels for StuPPOl, 
StuPP02, StuPP03 and StuPP04 (Figures 3D), and their 
similar low values for RPR, RPPO and RBR (Figure 4A, B 
and C). The three amiRPP0234A lines further clustered 
with line amiRPP023-5, because the four lines performed 
similarly in almost all of the variables except that the 
StuPPOl mRNA level in the amiRPP023-5 was moderately 
higher than in the amiRPP0234A lines (Figures 3C, D, and 
4). Three amiRPP01234 series lines (-2, -6 and -12) 
grouped closely because they had very similar expression 
trends in all of the 7 variables (Figures 3E and 4). 



StuPP04 gene expressions and had the least browning po- 
tential, were grouped within the lower quadrant with nega- 
tive factor scores in both PCI and PC2, and opposite the 
WT with positive factor scores in both PCI and PC2. 



Discussion 

PPOs are encoded by a gene family composed of mul- 
tiple highly conserved gene members in many plant 
species [2,24]. The differential temporal and spatial 
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Figure 5 Browning phenotypes. Potato tubers from each transgenic line and wild type were randomly selected, cut into approximately two 
equal sections and exposed to air/oxygen at room temperature (~25°C). Images were taken at 0, 24 and 48 hour post air exposure. Browning or 
blackening typically started from vascular ring region and advanced to the medullary tissue with the air-exposure time. 
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Table 4 Pearson correlation coefficients analysis results 



RPR RPPO RBR 



StuPPOl 


0.20 


0.19 


0.27 d 


SWPP02 


0.70 a 


0.78 a 


0.80 a 


SWPP03 


0.59 a 


0.69 a 


0.71 a 


StuPP04 


0.66 a 


0.69 a 


0.65 a 


RPR 




0.85 a 


0.89 a 


RPPO 






0.87 a 



Note: Pearson correlation coefficients (r 2 ) of the relationship of PPO gene 
expression level {StuPPOl, StuPP02, SWPP03 and SWPP04), relative PPO 
protein expression level (RPR), relative PPO enzymatic activity (RPPO) and 
relative browning potential (RBR) in potato tuber tissues, a, P < 0.0001; 
d, P < 0.05. 



expression patterns of PPO genes in potato, poplar and 
tomato indicate the functional diversities among the 
PPO gene members [1,33,34]. However, the lack of 
specific loss-of-function mutants of different PPO 
genes has impeded progress in understanding the di- 
versified gene functions in the family. Before the 
present report, several studies reported knockdown of 
PPO gene expression using siRNA strategies by trans- 
formation with hundreds of base pairs of the PPO gene 
fragments [4-10]. The strategies proved successful in 
suppression of overall PPO activities and reduction of 
PPO-mediated browning reactions, but the methods 
were not able to identify the roles of the individual 
PPO gene members. The heterogeneous siRNAs 
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Figure 6 Score scatter plot of the transgenic lines and WT 

according to principal component analysis (PCA). The bioplot 

represents the first and second principle components (PC1 and PC2, 

eigenvalue > 1.0) produced from the PCA of 15 observations (14 
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generated in the transformed plants did not target a 
specific PPO gene but affected others of the family as 
well because the PPO gene family shares high hom- 
ology in nucleotide sequences among its members 
[2,24]. Our data demonstrated that PPO gene isoforms 
can be suppressed individually or simultaneously using 
amiRNAs. In addition to significantly down-regulating 
the PPO genes of interest as predicted in our construct 
design, several exceptions were observed based on the 
qRT-PCR assay results (Figure 3). For example, an un- 
expected moderate to high reduction of the StuPP02 
and StuPP03 mRNA levels in line amiRPPOl-2, a 
StuPPOl gene knockdown transformant, and an un- 
foreseen moderate decrease of the StuPPOl and 
StuPP02 gene transcript copies in line amiRPP03-12, 
a StuPP03 gene knockdown mutant were detected al- 
though the analogue regions of the affected genes have 
6-8 bp mismatches to the targets of the designed 
amiRNAs (Figure 3A and B, Additional file 1: Part D). 
An 8-bp complementarity to the 5' seed region (ca- 
nonical 8mer site) [35] of the amiRPP03 in a StuPP04 
gene sequence region might be related to an unantici- 
pated moderate down- regulation of the non-targeted 
StuPP04 gene in Line amiRPP03-15 (Figure 3B, Add- 
itional file 1: Part D). Apart from strong down-regulations 
of both StuPP02 and StuPP03 genes in lines 
amiRPP023-5, -7 and -9, StuPP04 gene mRNA in these 
lines was at a moderate to low level, suggesting StuPP04 
gene expression was also targeted by amiRPP023 (Fig- 
ure 3C). The amiRPP023 was initially designed to target 
both StuPP02 and StuPP03 genes based on incomplete 
StuPP04 gene sequence data (GenBank ID: U22923). 
Later analysis of a StuPP04 gene related EST sequence 
(GenBank ID: BG592710) found a region with only one 
mismatch to the target of amiRPP023 (position 11, Add- 
itional file 1: Part D). A previous study reported mis- 
matches at position 11 were miRNA admissible targets 
but reduced miRNA-guided cleaving efficiency in vitro 
and phenotypic effects in vivo [36]. Although the possible 
target region of the StuPPOl gene has mismatches at two 
positions compared to the perfect target sequence of 
amiRPP0234A (positions 6 and 9 at the amiRPP0234A 
site, Additional file 1: Part D), the StuPPOl gene was 
moderately down-regulated in amiRPP0234A-4, -6 and 
-14 lines, along with strong down-regulation of the 
StuPP02, StuPP03 and StuPP04 genes (Figure 3D). It is 
worth noting that G:U wobble pairing was introduced for 
targeting StuPP03 gene in amiRPP0234A and StuPPOl 
gene in amiRPP01234 (Table 2) [37], and the corresponding 
target genes, StuPP03 gene in the amiRPP0234A series 
lines and StuPPOl gene in the amiRPPOl 234 series lines, 
were indeed down-regulated in both cases (Figure 3D and 
E). We failed to detect the designed amiRNAs from any of 
the transgenic plants transformed with the constructs for 
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amiRPP02 (transgenic amiRPP02 series lines) and 
amiRPP0234 (transgenic amiRPP0234 series lines) for 
unknown reasons (Table 3). Investigation of the principles 
and complexities between amiRNAs and target recogni- 
tion was beyond the scope of this study. Nonetheless, our 
data showed that amiRNAs are a potent regulator for 
modulating expression of potato PPO gene isoforms indi- 
vidually or in combination, and amiRNA targets are more 
predictable and detectable than the siRNAs. However, not 
all designed amiRNAs were expressed with similar effi- 
ciency in the transgenic plants and off-target effects also 
occurred occasionally. Appropriate detection methods are 
required for screening the desired transgenic traits in 
plants expressing amiRNAs. 

Individual or cumulative down-regulation of PPO genes 
did not usually cause up- regulation of the other PPO 
genes, except for transgenic line amiRPPOl-12 that 
showed a reduction of the targeted StuPPOl gene 
mRNA level by about 90% but also an increase of 
approximately 40% of the non-targeted StuPP02 gene 
transcript (Figure 3A). This implied that different 
PPO genes may be regulated independently reflecting 
their diversified functions in the potato plant. Sup- 
pression of PPO gene expressions in our transgenic 
plants resulted in a reduced level of total PPO protein, 
but the reductions varied depending on the knock- 
down targets. For example, expression of StuPPOl was 



suppressed by 90 to 98% in lines amiRPPOl -3 and 
amiRPPOl-12 (Figure 3A), and the total PPO protein 
level in these two transgenic lines was reduced by -20 to 
30% (Figure 4A). Coincidentally, the transgenic lines 
amiRPP023-5, -7 and -9, which simultaneously down- 
regulated the StuPP02, StuPP03 and StuPP04 gene 
(Figure 3C), had reduced total PPO protein level of about 
70% (Figure 4A). Total PPO protein level was reduced 
by 15% on average in line amiRPP03-15 that showed 
suppression of the StuPP03 gene (target) and StuPP04 
gene (off- target) (Figures 3B and 4A). These data indicated 
that StuPP02 protein probably contributed 55% or more 
to the total PPO protein levels in the non-transgenic po- 
tato tuber tissues. The failure to generate StuPP02-knock- 
down mutants made it difficult to evaluate the accuracy of 
the estimation. However, Pearson correlation coefficient 
indicated that 70% of the variation in total PPO protein 
level in potato tuber tissue was explained by StuPP02 
gene expression (r 2 = 0.70, Table 3). Analysis of PPO 
gene expression in the non-transgenic potato tubers by 
qRT-PCR found that the StuPP02 gene contributed 
67% to the total PPO gene transcript level, followed by 
StuPPOl, 28%, StuPP03, 4% and StuPP04, less than 
0.2% (Additional file 6: Table S2). This also indicated 
that StuPP02 is a major contributor to the PPO pro- 
tein content of the WT potato tubers. Noticeably, 
StuPPOl, StuPP02, StuPP03 and StuPP04 genes were 
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all down- regulated by 90% or more in transgenic 
amiRPP01234-2, -6 and -12 lines (Figure 3E), but the 
total PPO protein level remained between 20 and 30% 
in their tuber tissues. An explanation for this is that 
the PPO proteins that accumulate in the plastid prob- 
ably accumulated for a longer duration than the 
cytosol-produced PPO mRNAs [1,2]. 

A striking reduction of PPO-mediated browning has 
been described previously through suppression of overall 
PPO gene expression using siRNA strategies [4-6]. Here, 
we have shown that it requires that StuPPOl, StuPP02, 
StuPP03 and StuPP04 genes were all suppressed simul- 
taneously to achieve optimal inhibition of PPO activity 
and non-browning phenotype in transgenic potato tuber 
tissue (Figures 4 and 5), a point not well known or 
understood previously. It was of particular interest that a 
series of low browning phenotypes were generated 
though suppression of the highly conserved PPO gene 
isoforms when suppressed both individually and in 
combination using this amiRNA strategy (Figure 5). 
Analysis of these transgenic clones clearly demon- 
strated the positive correlations between PPO protein 
levels (RPR), PPO enzymatic activity (RPPO) and 
browning potential (RBR) (Figure 4, Table 4). For 
example, a -20-30% reduction at PPO protein level in 
the tubers of lines amiRPPOl-3 and -12, where the 
StuPPOl gene expression were down-regulated by over 
90%, showed a parallel decrease of both PPO activity 
and browning potential in the tuber tissues. Simultan- 
eous down-regulation of StuPPOl, StuPP02, StuPP03 
and StuPP04 genes in lines amiRPPO 1234-2, -6 and -12 
resulted in 70-80% loss of tuber PPO protein. As a result 
of this, PPO enzymatic activity was reduced by -90% and 
browning potential by 70-75%. Browning potential corre- 
lated better with the PPO protein level than with PPO en- 
zymatic activity. PPO enzymatic activity appeared to 
decrease relatively faster in our assays than browning reac- 
tion (Figure 4). The difference may have resulted from dif- 
ferent sample preparation methods used for the two 
assays. PPO enzymatic activity was more temporally sensi- 
tive after the enzymes were mechanically isolated from 
the plastids, while browning potential recorded the accu- 
mulation of pigments in the samples during the time 
period of the assay. In spite of the small discrepancy be- 
tween these two data sets, the trend was consistent among 
the PPO gene down- regulated transgenic plants; the lower 
the PPO protein level in the potato tuber tissues, the less 
the PPO activity and browning potential. 

Our survey of the recent S. tuberosum whole genome 
database (http://www.jgi.doe.gov) revealed 9 PPO-like 
gene models (Table 1, Additional file 1). StuPPOl to 
StuPPOS are clustered on chromosome 8 and only 
StuPP09 is located on chromosome 2. StuPPOl, StuPP02, 
StuPP03 and StuPP04 are allelic to previously characterized 



POTP1/P2, POT32, POT33 and POT72, respectively, in this 
report, based on their 95-99% identities at nucleotide level 
[1,12]. The deduced peptides also show 95-98% identical po- 
sitions between the corresponding alleles (Additional file 1). 
Surprisingly, StuPPOl is the only possible allele to the 
POTP1 and the POTP2 genes, and no duplication of the 
StuPPOl locus in the S. tuberosum genome was discovered. 
StuPPOS to StuPP09 are five novel PPO-like gene models 
predicted in the S. tuberosum genome (Table 1, Additional 
file 1). However, the low to no prevalence in the EST data- 
bases for the StuPPOS, StuPP06, StuPP07 and StuPPOS 
implies these gene models may be under-transcribed in S. 
tuberosum (Additional file 1: Part C). In addition, the im- 
portant tyrosinase domains (pfam00264) for PPO enzymatic 
activity are incomplete in the peptides of the deduced trans- 
lation of the StuPPOS and StuPP07 gene models, suggesting 
that these two may not be functional like other PPO iso- 
forms (Additional file 1: Part B) [23,24]. StuPP09 is probably 
a PPO-like gene involved in disease defense and cell rescue 
because the ESTs that support for this gene model were 
dominantly discovered in the tissues of in vitro cultured po- 
tato callus, abiotic stress treated leaf and root or pathogen- 
infected leaf and tuber, and several of the ESTs expressed at 
significantly higher levels in the pathogen-infected potato 
tissues (Additional file 1: Part C) [25-28]. Taken together, 
the number of members of the potato PPO gene family are 
possibly larger than previously reported, but StuPPOl 
(POTP1/P2), StuPP02 (POT32), StuPP03 (POT33) and 
StuPP04 (POT72) are probably the major developmental^ 
regulated PPO genes in S. tuberosum and they are subjected 
to the targets of the amiRNAs expressed in the transgenic 
lines. However, extensive experimental studies are required 
to investigate the existence and function of the potential 
StuPPOS to StuPP09 genes. 

Conclusions 

PPO-mediated browning damage is one of the main 
causes of quality loss in fresh fruit and processed food. It 
is of great importance and interest to produce crop var- 
ieties with low PPO activity for the food industry. Our 
results have shown that StuPPOl (POTP1/P2), StuPP02 
(POT32), StuPP03 (POT33) and StuPP04 (POT72) genes 
were the major contributors to the total PPO protein con- 
tent but the effect of the individuals was not equal and 
that PPO activity in the tuber tissues paralleled the protein 
content data. Suppression of expression of one or a few 
PPO genes did not cause overexpression of the others, but 
the greatest reduction of PPO activity and the most 
complete non-browning phenotypes required simultan- 
eous suppression of the expression of StuPPOl, StuPP02, 
StuPP03 and StuPP04 genes. Our demonstration that 
PPO gene expression in potato can be suppressed by 
introduction of 21-nt small RNA regulators provides a 
new strategy for developing low- or non-browning crops. 
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It also suggests that amiRNAs can be applied to silence 
closely related members of multi-gene family for func- 
tional genomics study in non-model plants. In addition, a 
series of PPO gene knockdown mutants predictably gener- 
ated provide us important resources for future investiga- 
tions of the role of PPO genes in functions such as 
plant development, stress and insect and fungal defense 
response. 

Methods 

Survey of PPO-like genes in Solanum tuberosum 

Nucleotide sequence information of the identified potato 
PPO genes, including POTP1 (GenBank ID: M95196), 
POTP2 (M95197), POT32 (U22921), POT33 (U22922), 
and POT72 (U22923) were obtained from the GenBank. 
The gene sequences were used as queries to search for 
the PPO-like gene sequences from a recent S. tuberosum 
whole genome database in the US Joint Genome Insti- 
tute (http://www.jgi.doe.gov) using the blastn with the 
default parameters. The BLAST hits were manually 
checked and the PPO-like gene sequences were retrieved 
and further analyzed using Vector NTI Advance 11 soft- 
ware (Life Technology, USA), NCBI BLAST platform 
(blastn, blastp, nr/nt, ESTs, etc.), and Simple Modular 
Architecture Research Tool (SMART, http://smart.embl- 
heidelberge.de/). 

Construction of amiRNA vectors 

The amiRNAs for targeting the characterized potato PPO 
genes, including StuPPOl (previously named P0TP1/P2), 
StuPP02 (previously, POT32), StuPP03 (previously, POT33) 
and StuPP04 (previously, POT72), individually or cumula- 
tively were designed in an Arabidopsis miR-168a gene back- 
bone using an in-house amiRNA designer and constructed 
as previously described [30]. All constructs were sequenced 
to confirm the intended construction and the designs. De- 
scriptions of the amiRNAs, their target sequences and genes 
are listed in Table 2. 

Plant genetic transformation 

Transformation of amiRNA constructs with potato culti- 
var ALT 1762 was performed as previously described 
[30]. Transgenic plants with each of the amiRNA designs 
were assigned a series of names and numbers as listed in 
Table 3. 

Analysis of PPO gene expression by real-time quantitative 
reverse transcription PCR (qRT-PCR) 

Potato tissues were disrupted in liquid nitrogen using 
mortar and pestle. Total RNA was extracted using a 
Spectrum plant total RNA kit (Sigma- Aldrich, USA) ac- 
cording to the manufacturer s instruction. The purified 
total RNA samples were treated with Ambion TURBO 
DNase (Life Technologies, USA) to eliminate potential 



contamination from genomic DNA, and subsequently 
re-purified using RNeasy mini kit (Qiagen, Germany). 
RNA concentration and purity were determined using a 
ND-1000 spectrophotometer (NanoDrop, Wilmington, 
DE, USA). First strand cDNA was synthesized from 1 ug 
of the re-purified total RNA using 1 ul of RT primer 
mix (provided with the QuantiTect reverse transcrip- 
tion kit, Cat. No. 205310, Qiagen) and 1 ul of Super- 
Script III reverse transcriptase (200 U/ul) in a 20-ul 
reaction according to the manufacturers instruction 
(Life technologies). 

Two housekeeping genes, translation elongation factor 
1 -alpha (efla) and cyclophilin were selected as refer- 
ences for qRT-PCR data analysis. The primer pairs for 
efla and cyclophilin gene amplicons were the same as 
the designs by Nicot et al. [38] (Additional file 7: Table 
S3). The expression stability of the two reference genes 
in potato tuber tissues was validated by geNorm algo- 
rithm [39] (Additional file 8: Figure S4). More than 30 
primers for amplions of PPO genes were designed using 
the Primer 3 (version 0.4.0) [40] based on available 
gene-specific nucleotide sequence information retrieved 
from the GenBank databases. Four primer pairs, one for 
each of the PPO genes, namely StuPPOl, StuPP02, 
StuPP03 and StuPP04 were eventually selected for 
qRT-PCR analysis in this study. The PCR amplification 
specificities using these primer pairs were examined by 
agarose gel electrophoresis, cloning and sequencing, as 
well as post- PCR melting curve analysis (Additional file 
7: Table S3). All primers used in this research were syn- 
thesized by the Integrated DNA Technologies, Inc. 
(Iowa, USA). 

All qRT-PCRs were carried out in triplicate in 96- 
well microplates and performed in a CFX96 real-time 
PCR detection system (Bio-Rad, USA). Each reaction 
volume was 20 ul, comprised of 10 ul of Eva green 
master mix (Applied Biological Materials Inc., Canada), 
0.5 uM of forward and reverse primers (primer pair 
mix, 5 ul) and 5 ul of appropriately diluted DNA or 
cDNA template. No-template control (NTC) reactions 
were included in each plate to monitor potential forma- 
tion of primer-dimers. The thermal cycling was pro- 
gramed as follows: one initial denaturation cycle at 95° 
C for 10 min; 40 cycles at 95°C for 30 s and 60°C for 
30 s. Fluorescence signal was measured at the end of 
each annealing and extension step at 60°C. At the end 
of the qRT-PCR run, a melting curve analysis with a 
temperature gradient of 0.1°C/s from 65 to 95°C was 
performed to ensure that only single products were 
produced. Relative expression levels of the PPO gene 
expressions were calculated based on a method de- 
scribed in Pfaffl [41], Data normalization was per- 
formed using the gene expression values of the efla 
and cyclophilin in the samples. 
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Detection of amiRNAs by reverse transcription PCR (RT-PCR) 

Detections of mature amiRNAs expressed in the trans- 
genic plants by RT-PCR were conducted as previously 
described [30]. The forward primers used in the PCR 
were identical to the corresponding amiRNA sequences 
except that they were synthesized as oligodeoxynucleo- 
tides (Additional file 9: Table S4). 

Detection of cleaved mRNA of amiRNA target 

The FirstChoice RLM-RACE kit (Life Technologies) was 
explored to detect the cleaved mRNA product of the 
amiRNA target. The kit is designed to amplify cDNA 
only from full-length, capped mRNA by removing free 
5'-phophates from molecules such as fragmented mRNA 
using Calf-intestinal alkaline phosphatase (CIAP) before 
processing to remove the cap structure from full-length 
mRNA using Tobacco acid pyrophosphatase (TAP). In 
our application, the CIAP and TAP treatment steps were 
omitted to prevent full-length, capped-mRNA from ligat- 
ing to the 5 '-RACE adapter (Figure 2A). 

Poly(A) + RNA was purified using the Dynabead mRNA 
Purification kit (Life Technologies). Enriched poly(A) + 
RNA was ligated to the 5 '-RACE adapter provided in the 
FirstChoice kit without treatment with CIAP and TAP for 
exclusion of all 5 '-capped mRNA. The conditions for 1st- 
strand cDNA synthesis, the followed-up PCR and nested 
PCR were carried out as instructed by the manufacturer. 
Besides that, the following primer strategy was developed 
for detecting possible fragmented mRNA from StuPPOl, 
StuPP02, StuPP03 and StuPP04 genes (Figure 2A). A re- 
verse primer designed from a conserved sequence region 
of all four PPO genes was used to synthesize the 1st- 
strand cDNA. Two other reverse primers designed from 
another two conserved regions of the four PPO genes 
were used with the common 5 '-RACE primer to conduct 
the lst-round PCR and the follow-up nested-PCR-1. 
Using products from the nested-PCR-1 as the template, 
nested-PCR-2 reactions were performed with the com- 
mon 5 '-RACE primer and four gene-specific reverse 
primers, respectively for detecting the bio-origination of 
the cleaved mRNA products. Distinct bands of appropri- 
ate size of PCR fragments were separated by electrophor- 
esis, cloned and sequenced. The oligonucleotides used for 
detection of cleaved mRNA are listed in Additional file 10: 
Table S5. 

Analysis of PPO protein levels in potato tuber tissues 

A protein dot blot assay was used for semi- quantifying 
PPO protein levels in potato tubers. Crude protein sam- 
ples were prepared by powdering 0.2 g tuber tissue in li- 
quid nitrogen followed by adding 1.6 ml of extraction 
buffer (2% SDS, 0.1% (3-mercaptoethanol in 63 mM 
Tris-HCl buffer, pH 6.8). Samples were boiled for 
10 min, followed by centrifugation at 20,000 g for 2 min 



to remove insoluble material. The crude protein concen- 
trations were determined using the BCA Protein Assay 
Kit (Thermo Scientific Pierce, USA) and all samples 
were then adjusted to 100 mg/ml protein stock. For dot 
blot assay, 1:400 diluted samples from the protein stocks 
were blotted onto polyvinylidene difluoride membrane 
using Easy-Titer ELIFA System 77000 (Pierce) and 
probed with a PPO specific antibody (Ab). The PPO 
Ab was raised in rabbit against a synthetic peptide 
(KDWLNSEFFFYDE) corresponding to a conserved re- 
gion in the deduced StuPPOl, StuPP02, StuPP03 and 
StuPP04 protein sequences by Applied Biological Ma- 
terials Inc. (Canada). Two-fold serial dilutions of the 
synthetic peptide (KDWLNSEFFFYDE) from 50 ng to 
0.78125 ng were included in each blot for internal con- 
trol and making protein standards. The PPO Ab probed 
blot was incubated with HRP- conjugated anti-rabbit IgG 
and detected with the Enhanced Chemiluminescence De- 
tection System (GE-Amersham, USA). PPO Protein con- 
centrations were determined using Aida Image Analyzer 
version 2.00. Relative protein expression level, i.e. PPO 
proteins in the transgenic plants to that in wild type was 
calculated for presentation in this study. 

PPO enzymatic activity assay 

PPO enzymatic activity was measured using 4-methylcatechol 
(4-MC) as substrate. A 40 mg sample with a dimension of 
2 mm x 5 mm x 5 mm was collected from the potato 
tuber tissue next to the skin using a fruit peeler and a 
single-hole paper punch, and frozen in liquid nitrogen im- 
mediately. The sample was homogenized with cold PPO 
extraction buffer (100 mM sodium phosphate buffer 
pH 6.0, 2% TX-100, 2% PVPP) using a Fast Prep FP120A- 
115 Homogenizer (Thermo Electron Corp, USA). A ten- 
fold diluted extract (20 ul) was mixed with 200 ul of the 
assay mixture (50 mM sodium phosphate buffer pH 6.0, 
0.1% SDS and 15 mM 4-MC) in 96-well Microplates. A 
SpectraMax M2 Microplate Reader (Molecular Devices, 
USA) was used to measure the absorbance increase at 
400 nm (A 40 o n m) every 5 seconds for 1 min at 25°C. One 
unit of PPO enzymatic activity was defined as the amount 
of enzyme necessary to change A 400nm in 0.001 /min at 
25°C. The total protein concentration of the extract 
was detected using the BCA Protein Assay Kit (Pierce) 
and measured in the same SpectraMax M2 Microplate 
Reader. The enzymatic activity was calculated as U/mg 
of total protein. 

Browning potential assay 

Browning potential of phenolic compounds in potato 
tuber tissues was measured using a modification of the 
method described previously by [32]. Each sample tissue 
was a collection of two randomly selected potato tubers 
harvested from the same transgenic line or wild type. 
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Triple biological replicates (representing a total of 6 tu- 
bers per line) were set up for each treatment. The skins 
were peeled off before the tubers were chopped into 
small pieces manually. A 3-gram subsample of each 
pooled two-tuber replicate was weighed out and homoge- 
nized with 15 ml cold sodium phosphate buffer (50 mM, 
pH 6.0) using an Ultra-Turrax tissue homogenizer 
(Takmar, Cincinnati, Ohio) at high speed for 1 min. 
The homogenate was allowed to oxidize at room 
temperature for 1 h before it was vortexed and 1 ml 
aliquots were centrifuged at 12,000 g for 10 min. 
300 \i\ of the supernatant was loaded in a 96-well mi- 
croplate for measuring the absorbance at 475 nm 
using a SpectraMax M2 Microplate Reader (Molecular 
Devices) and the A 475nm value was used as an indica- 
tion of the browning potential. 

Statistics analysis 

Analysis of variance with the general linear model proce- 
dures, means comparison by Duncans multiple range 
test and the Pearsons correlation analysis were per- 
formed according to Statistical Analysis System, SAS 9.1 
for windows (SAS Institute, Cary, NC). Principal compo- 
nent analysis (PCA) was conducted with XLSTAT 2012 
for Windows (Addinsoft, NY). Hierarchical cluster ana- 
lysis (HCA) was carried out with SYSTAT 12 version 
12.02 for Windows (Addinsoft, NY) using Wards link- 
age hierarchical clustering with the Euclidean distance. 

Additional files 



Additional file 10: Table S5. List of reverse primers used in detection 
of the cleaved mRNA of amiRNA target in transgenic lines amiRPP01-12 by 
5 -RACE PCR. 
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